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We theoretically study thermally-activated phase slips in superfluid spin transport in easy-plane 
magnetic wires within the stochastic Landau-Lifshitz-Gilbert phenomenology, which runs parallel 
to the Langer-Ambegaokar-McCumber-Halperin theory for thermal resistances in superconducting 
wires. To that end, we start by obtaining the exact solutions for free-energy minima and saddle 
points. We provide an analytical expression for the phase-slip rate in the zero spin-current limit, 
which involves detailed analysis of spin fluctuations at extrema of the free energy. An experimental 
setup for a magnetoeletric circuit is proposed, in which thermal phase slips can be inferred by 
measuring nonlocal magnetoresistance. 

PACS numbers: 75.76.-hj, 74.20.-z, 75.78.-n, 75.10.Hk 


Introduction .—A wire can carry an electrical current 
without dissipation under favorable conditions in the su¬ 
perconducting state, which is characterized by a complex- 
valued function of position, ^(r), referred to as the su¬ 
perconducting order parameter describing a condensate 
of constituent particles [1]. When density fluctuations of 
the condensate are energetically suppressed and thus the 
magnitude of the order parameter is constant, an elec¬ 
trical supercurrent is proportional to the gradient of the 
phase of the order parameter. In some circumstances, 
e.g., for thin wires or in the presence of strong magnetic 
fields, finite resistances arise [2], of which understand¬ 
ing has required both theoretical and experimental ef¬ 
forts to be made over the last decades. In particular, 
the theory for intrinsic thermal resistances in thin su¬ 
perconducting wires has been pioneered by Little [3] who 
suggested that dissipation occurs via thermally-activated 
phase slips (TAPS), jumps of the phase winding along the 
wire by 27r, that necessarily accompanies quenching of the 
order parameter at some point in the wire. The quan¬ 
titative theory for TAPS has been developed by Langer 
and Ambegaokar [4] and further elaborated by McCum- 
ber and Halperin [5], and has, therefore, been given the 
acronym following their initials: LAMH theory. 

No superconductivity has been observed at room tem¬ 
perature, challenging its practical utilization, whereas 
magnetism—another phenomenon resulting from spon¬ 
taneous ordering—is ubiquitous in nature even at ele¬ 
vated temperatures. Being integrated with the infor¬ 
mation processing technology, it has spawned the field 
of spintronics [6]. A spin analog of an electrical super¬ 
current, superfiuid spin transport, has been proposed in 
magnets with easy-plane anisotropy, where the direction 
of the magnetic order parameter within the easy plane 
plays a role of the phase of the superfluid order param¬ 
eter [7-10]. Here, dissipationless spin current (polarized 
out of the easy plane) is sustained by a planar spiral¬ 
ing texture of the magnetic order. The absence of strict 
conservation laws for spin, e.g., due to Gilbert damp¬ 
ing, rules out faithful analogy to electrical supercurrent. 


which requires us to consider superfiuid spin transport as 
being distinct from conventional charge superfiuid. 

In this Letter, we theoretically study TAPS in super- 
fluid spin transport in easy-plane magnetic wires within 
the Landau-Lifshitz-Gilbert (LLG) phenomenology. In 
equilibrium, the magnetic order is kept in the easy plane 
and, thus, can be characterized by its winding number, 
the total azimuthal-angle (phase) change along the wire. 
At a finite temperature, the winding number can increase 
or decrease due to thermal spin fluctuations via events 
that can be identified as TAPS. The most probable path 
for the dynamics of the magnetic order parameter dur¬ 
ing TAPS traverses the saddle point of the free energy, 
where a few localized spins develop significant out-of- 
easy-plane components. We obtain the exact solution 
for these saddle points by solving the time-independent 
Landau-Lifshitz equation. We also provide an analytical 
expression for the rate of TAPS in the zero spin-current 
limit, which involves detailed analysis of spin fluctuations 
at extrema of the free energy. To observe TAPS in mag¬ 
netic wires, we adopt a magnetoeletric circuit proposed 
in Ref. [11], in which detection of nonlocal magnetoresis¬ 
tance can yield signatures of TAPS. 

Main results .—We consider a thin easy-x^-plane mag¬ 
netic wire with free energy F[n] = f dV{A\'Vn\‘^ + 
Kn‘l)l2, where positive constants A and K parametrize 
the stiffness of the order parameter and the easy-plane 
anisotropy, respectively. Here, the unit vector n{x) is the 
direction of the order parameter: the local spin angular- 
momentum density for ferromagnets and the local Neel 
order for antiferromagnets. When the wire is narrow 
compared to the magnetic coherence length ^ = ^/A/AT, 
variations of the order parameter across the wire (of cross 
section S) can be neglected, which allows us to treat 
the order parameter n, at a given time, as a function 
of position x along the wire [12]. It is convenient to 
parametrize n in spherical coordinates, 0 and 0, defined 
by n = (sin6>cos0, sin6>sin0, cos6>), with a rescaled free 
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FIG. 1. (color online) (a)-(c) Extrema of the free energy / 
[Eq. (1)] for an easy-plane magnetic wire with periodic 
boundary conditions. The metastable state (a) that winds 
in the easy plane once will decay to the ground state (c) via 
TAPS passing over the saddle point (b), where a few spins 
localized within the magnetic coherence length develop sig- 
nihcant out-of-easy-plane components, (d) A plot of the free 
energy / as a function of spin current k (metastable states) 
and k (saddle points) for the wire length / = 48. A solid line 
is a guide to the eye. A dashed line shows the free energy of 
the metastable states for an inhnitely long wire: / = k‘^l/2. 
Points corresponding to the three conhgurations, (a)-(c), are 
denoted accordingly. A dotted line illustrates transitions be¬ 
tween nearby metastable states (a) and (c). 


energy: 

nlt2 

f[0^(j)\= / dx[0'^0(j)'^ ^ 0]/2 (1) 

J-112 

measured in units of Fq = S^K (which is the maximum 
anisotropy energy that can be stored within the coher¬ 
ence length), where the position variable x is measured 
in units of ^ and runs over the wire length I = L/^. 

Extrema of the free energy are solutions of the time- 
independent Landau-Lifshitz equation: 

Sf /60 = —0" + sin 0 cos 0(j)'^ — sin 0 cos ^ = 0 , (2a) 

8f/5(j)=—{^\v?0(j)'y = ^. (2b) 

The second equation is the consequence of the invariance 
of the free energy under spin rotations about the ^ axis. 
For static configurations, the associated conservation law 
describes spatial independence of the 2 :-component of the 
spin current, Ig = — AF sin^ so that the dimension¬ 
less constant parameter k = —Isi/AS can be used to 


index solutions of Eqs. (2). There are two types of solu¬ 
tions of interest to us. The first is a local minimum of 
the free energy: 


0{x)=tt/2, (j){x) = (t)o ^ kx (|/c|<l), (3) 

with 00 an arbitrary reference angle. There is a crit¬ 
ical current, \k\ = 1, for stable superfiuid spin trans¬ 
port according to the Landau criterion [7, 8, 13], above 
which spin fluctuations destabilize superfluidity. When 
the wire is long enough I ^ 1 (which is assumed hence¬ 
forth), actual boundary conditions at the ends of the 
wire are not important. Imposing periodic boundary 
conditions on the order parameter, n{x = —1/2) = 
n{x = //2), quantizes the total azimuthal-angle change: 
A0 = 0(//2) — 0(—//2) = 27r/i, in terms of integer ji. 
The allowed values of k are thus k^ = 27r/i//. Figures 1(a) 
and 1(c) show the free-energy minima with winding num¬ 
bers /i = 1 and /i = 0, respectively. 

At zero temperature, thermal spin fluctuations are 
frozen out. Persistent spin current in a closed magnetic 
ring, therefore, can be sustained indefinitely, when dis¬ 
regarding quantum spin fluctuations [14]. Finite tem¬ 
perature, however, agitates spins and opens transition 
channels between the metastable states carrying differ¬ 
ent spin current [see a dotted line in Fig. 1(d)]. The 
total azimuthal-angle change A0 = 27r/i is quantized and 
well defined provided that the order parameter n avoids 
the poles, \nz\ = 1, where the azimuthal angle 0 is am¬ 
biguous. In continuous transitions between two minima 
with different winding numbers, /i 7 ^ /i', the order param¬ 
eter must hit one of the poles; this is analogous to the 
vanishing of the superconducting order parameter during 
TAPS [3]. Supposing T ^ Fq, the transitions between 
metastable states are rare, which we assume throughout. 

The most probable path of the order parameter during 
the transition between two metastable states will pass 
over the intervening saddle point of the free energy [2, 5], 
which is the second kind of solution of Eq. (2) that we 
obtain with spatially varying 0{x) [15]: 


0{x) = cos ^ \/l — sech(\/1 — k‘^x) 


(4a) 


0(x) = (po kx tan 


-1 


tanh( \/l — k^x) 


(4b) 


indexed by spin current k, and any spatial translation 
thereof. This exact saddle-point solution constitutes our 
first main result. Periodic boundary conditions on n dis¬ 
cretize allowed values of k: A0 = k^l + 2tan“^[(l — 
= 27r/i, where the quantities exponentially 
small for large I are ignored here and hereafter. Fig¬ 
ure 1(b) depicts the saddle-point solution with /i = 1, 
which mediates the transition between two minima with 
/i = 1 and /i = 0. The spin currents of the metastable 
states and the saddle-point solutions interlace: k^-i < 














3 


kfx < (for positive /i), meaning that there always exists 
the unique saddle point between two nearest metastable 
states. See Fig. 1(d) for an illustration. 

The rate of transitions, respectively increasing or de¬ 
creasing spin-current magnitude, may be written in the 
form 

r± = , (5) 

where temperature is measured in energy units so that 
= 1. Here, AF± = Fq • A/± is the free-energy barrier 
to reach the intermediate saddle point, and Vi is the pref¬ 
actor that depend on details of spin fluctuations around 
the extrema [16]. Specifically, for the transitions between 
the two metastable states [Eq. (3)] with and k^-i via 
the saddle point [Eq. (4)] with k = k^ > 0, the free- 
energy barriers can be directly obtained by evaluating 
the differences in the free energy / [Eq. (1)]: 

Af_(k) = 2 \/1 — k‘^ — 2^tan“^[\/1 — k‘^/k] , (6a) 

A/+(^) = Af-(k) + 27Tk . (6b) 


Since A/_ < A/+, fluctuations tend, on average, to re¬ 
duce the spin-current magnitude and thus give rise to 
equilibriation. In the limit of zero current, ^ ^ 0, the 
free-energy barrier is AF = 2Fo = 2S^K, which roughly 
represents the energy cost due to the out-of-easy-plane 
component of the order parameter in the phase slip re¬ 
gion localized within the magnetic coherence length 
Our second main result, which is derived in the sup¬ 
plemental material [17], is the analytical expression of 
the prefactor Q for ferromagnets in the zero spin-current 
limit: 


n{T) 


1 aK L rKF 
7rv^(l + a2)siV^’ 


( 7 ) 


which is analogous to the result for the superconduct¬ 
ing wire in the LAMH theory [5], where a is the Gilbert 
damping constant and s is the local spin angular mo¬ 
mentum density. Here, aK/{l + a‘^)s is the inverse of 
the relaxation time for the perturbed uniform easy-plane 
ferromagnet to return to the equilibrium state; L/^ rep¬ 
resents the number of possible independent phase-slip lo¬ 
cations; AF/T stems from the breaking of the transla¬ 
tional invariance of the system by the saddle point [18]. 
The prefactor for antiferromagnets on bipartite lattice 
can be obtained by replacing aK/{l + (a^)s with K/as 
for overdamped dynamics [19], where s is the local spin 
angular-momentum density per each sublattice. 

Decay of persistent spin current .—The persistent spin 
current in a closed ring will decay via TAPS at a fi¬ 
nite temperature. Erom Eq. (5), the winding number 
fi = A0/27r, which characterizes metastable states, de¬ 
cays with the rate 

r+ - r_ = -47r2(^Fo/Lr)^^(T)e-2^“/^//, (8a) 

= -k{T)ii (8b) 
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FIG. 2. (color online) (a) Schematics of an experimental setup 
for detecting TAPS, in which two identical metals, parallel in 
the electric circuit, are connected by a magnetic insulator sup¬ 
porting superfluid spin transport, (b) Schematics illustrating 
the origin of an electromotive force in the metals. TAPS un¬ 
wind the equilibrium spiraling structure (at t = ti), resulting 
in the uniform state (at t — t 2 ). As the magnet returns to the 
equilibrium spiraling structure, the magnetization at the left 
(right) interface rotates counterclockwise (clockwise), which 
in turn induces a detectable electromotive force in the metals. 


to linear order in the winding number fi [20]. The 
spatially-averaged spin current A = 2'K[iASIL decays 
with the rate n{T)Is. Note that n{T) is independent of 
the length of the wire since fl(T) oc L. 

The dissipation of the spin current dictates the pres¬ 
ence of the effective random force on the spin current to 
meet the fluctuation-dissipation theorem [21]. The resul¬ 
tant stochastic dynamics of the spin current is described 
by 


is{t) = -K{T)Is{t)Fn{t), (9) 

where the white-noise Langevin term u(t) with the cor¬ 
relator {n{t)n{t')) = 2{AS/L)K,{T)TS{t — t') is intro¬ 
duced to yield the thermal variance of the spin current, 
(/^) = {AS/L)T, which we obtain from the thermal ex¬ 
pectation value of the free energy. 

Discussion .—TAPS in superfluid spin transport can be 
detected in an experimental setup proposed in Ref. [11], 
in which two identical metals connected parallel in the 
external electric circuit are linked by a thin easy-plane 
magnetic insulating wire (see Fig. 2). In the presence of 
spin-orbit coupling at metal| magnet interface, current in 
the metal gives rise to a torque in the magnet, and, as an 
Onsager reciprocal effect, dynamics of magnetic moments 
induces an electromotive force in the metal [22]. 

At zero temperature, this configuration supports static 
spiraling structure of the magnetization [11], with the left 
metal injecting and the right metal draining spin current 
[see Fig. 2(a)]. The associated electromotive force is ab¬ 
sent, and the effective resistivity of the circuit, therefore, 
is not affected by spin superfluid. At a finite temperature, 
however, TAPS unwind the spiraling structure stochas¬ 
tically with the net rate of /i:(T)/i(t), where /i(t) is the 
winding number at fixed time t. As the magnetization 
rewinds to the equilibrium spiraling structure, the mag¬ 
netic moment at the left (right) interface rotates coun¬ 
terclockwise (clockwise), which induces an electromotive 
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force in the adjacent metals [23]. In the steady state, 
the time-averaged precession of the magnetization at the 
interfaces are (pi{t) = —(j)r{t) = —7rK{T)jii{t). Induced 
electromotive force reduces the effective resistivity of the 
circuit (following the derivation of Ref. [11]): p ^ pm 
with pm = where p is the re¬ 

sistivity of the metal, is the effective interfacial spin¬ 
mixing conductance (parametrizing the Gilbert-damping 
enhancement at the interfaces due to spin pumping), and 
'd is related to the effective interfacial spin Hall angle 0 
via'd = (/i/2et) tan 0, with t being the thickness of the 
metals in the x direction [11]. Observation of the charac¬ 
teristic dependence of the effective resistivity of the cir¬ 
cuit on the length of the wire (algebraic) or temperature 
(exponential) would provide experimental signatures of 
TAPS. 

Spins are treated classically in our theory for TAPS. 
Quantum aspect of spins would become important in the 
low-temperature regime, where quantum phase slips may 
become a dominant source of dissipation of spin super¬ 
current. It is possible to pursue a study of quantum phase 
slips with the aid of semiclassical quantization of spins, 
which has been used to understand macroscopic quantum 
tunneling of magnetic domain walls [14]. Recalling about 
the superconductor-to-insulator quantum phase transi¬ 
tion (as a function of wire cross section) [2, 24], a natural 
question arises on the role of the spin Berry phase in the 
analogous physics of the magnetic quantum phase slips. 
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by the U.S. Department of Energy, Office of Basic En¬ 
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(an SRC STARnet center sponsored by MARCO and 
DARPA). 
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In this supplemental material, we derive the prefac¬ 
tor CL(T) in the phase-slip rates in the limit of zero spin 
current. 

Phase-slip rates. —To derive the prefactor ^2(T), it is 
convenient to impose antiperiodic boundary conditions 
on the order parameter: n{x = —112) = —n{x = 1/2), 
with which a unique saddle point solution is present at 
k = 0. There are two degenerate ground states with 
the total azimuthal angle changes A(f) = Ttt. The sad¬ 
dle point passed during the transition between these two 
states is analogous to Eq. (4) with ^ = 0: 

0 = cos“^(sechx), 0 = 00 + {^/2)H{x ), (SI) 

where H{x) is the Heaviside step function. 

For the probability per unit time of a single TAPS 
event, adapting the derivation for superconducting wires 
in Ref. [SI] to our problem, we obtain the following ex¬ 
pression for general spin current k: 

r{k, T) = n{k, , (S2) 

with the prefactor 



FIG. SI. (color online) Extrema of the free energy [Eq. (1)] 
with antiperiodic boundary conditions. Transitions between 
the two metastable states, (a) and (c), occur via TAPS passing 
over the saddle point (b), where a few spins in the phase slip 
region localized within the magnetic coherence length develop 
significant out-of-easy-plane components. 


n{k,T) 


j{k)A{k) I Fq I det'PL 
^ V \ det'n\ ' 


(S3) 


Similar expression has been derived in Refs. [S 2 ] to obtain 
rates of magnetization reversal in ferromagnetic wires 
with easy-axis anisotropy, where transitions between two 
ground states with uniform magnetization are driven by 
soliton-antisoliton nucleations. Here, j{k) is the escape 
rate from the involved saddle point, according to the lin¬ 
earized equation of motion at the saddle point solved 
by 6n{x,t) = exp{^t)6n{x) with 7 = SPy [S3]; A{k) 
stems from the presence of the zero mode associated with 
the translation symmetry breaking by the saddle point; 
PL and PL (which are discussed further below) are the 
differential operators accounting for contributions from 
fluctuations Sn to the free energy / at the metastable 
states and saddle-point solutions, respectively; det^ O is 
the product of nonzero eigenvalues of the operator O. 

Fluetuations. —Let us investigate spin fluctuations at 
extrema of the free energy. To describe small devia¬ 
tions (5n from the extrema ne, it is convenient to setup 
the local ortho normal frame defined by 63 = ng, ei = 
dwe/dO, 02 = (sin^)“^9ne/90, in which the deviation is 
parametrized by two numbers u and v: n = uei -\- ve 2 + 
{1 — — i;^)^/^e3. 


The contribution from spin fluctuations, T(x) = 
[u{x),v{x)], to the free energy / at the metastable state 
with the spin current k is f dx T(x)^HT(x)/ 2 , where 


f—cP/dx‘^ 1 — k‘^ 0 \ 

0 —cP/dx‘^J 


(S4) 


The corresponding eigenvalue problems for u{x) and v{x) 
are decoupled: 

{—d^/dx‘^ + 1 — k^)un{x) = enUn{x ), (S5a) 

{-(T/dx^)Vn{x) = e'nVn{x). (S5b) 


Note that above the critical current, \k\ = 1, the stable 
spin current cannot be sustained due to unrestricted spin 
fluctuations. 

Similarly, the contribution from spin fluctuations to 
the free energy at the saddle point with k = 0 is 
f dx ^(x)^PL^(x)/2, where 

y _ f—^/dx^ + 1 — 2 sech^ x 0 

0 —— 2 sech^ X 

(S 6 ) 

The eigenvalue problem for u{x) and v(x) also decouples 
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[S4]: 

{—S/dx^ + 1 — 2 sech^ x)un{x) = enUn{x ), (S7a) 

{—d^/dx‘^ — 2 sech^ x)vn(x) = e^Vn{x). (S7b) 

Two eigenvalue problems, Eqs. (S5) and (S7), at /c = 0 
are related by the supersymmetry [S5] and are exactly 
solvable [S 6 ]. Two properly-normalized localized solu¬ 
tions of Eq. (S7) are of particular interest to us: 

e'l = — 1 , vi{x) = 2 “^/^ sechx , (S 8 a) 

ei = 0 , 'Ui(x) = 2 “^/^ sechx . (S 8 b) 

The negative-eigenvalue solution vi{x) represents the di¬ 
rection of deviation from the saddle point that de¬ 
creases the free energy the most rapidly, which shall be 
used to derive 7 in Eq. (S3). The zero-eigenvalue solution 
ui (x) is the Goldstone mode associated with the freedom 
of choice for a phase-slip location, which shall be used to 
derive A in Eq. (S3). 

Derivation of the prefaetor .—Let us now derive the es¬ 
cape rate from the saddle point, 7 , in n{T) [Eq. (S3)]. 
We will show details for the case of ferromagnets. The 
dynamics of n is described by the LLG equation, s{l -\- 
(anx)h = n X h, where h = —dF/dn is the effec¬ 
tive field conjugate to n. Writing ['u(x,t),f;(x,t)] = 
ex.p{jt)[u{x),v{x)], the linearized LLG equation at the 
saddle point for u{x) and v{x) is given by 

s^[u{x) — av{x)] = K{Hu — l)v{x ), (S9a) 

s^[v{x) -h au{x)] = —KHuu{x ), (S9b) 

with the differential operator 1-Lu = —d?ldx‘^ + 1 — 
2sech^x. It is easy to verify that the equations are sat¬ 
isfied by the escape rate 7 = aK/(l (y‘^)s in conjunc¬ 
tion with the eigenfunction [u{x),v{x)] = sech(x) [ 1 , —a] 
whose special profile was taken from Eq. (S 8 ). Applying 
the same method to antiferromagnets on bipartite lat¬ 
tice, whose dynamics is described by the antiferromag¬ 
netic LLG equation [S5, S7], we obtain 7 = K/as for 
over damped dynamics. 

A in Eq. (S3) can be derived following Ref. [SI]. In¬ 
finitesimal spatial translations of the saddle-point solu¬ 
tion [Eq. (SI)] can be represented by d0{x) = dx sech(x). 
Equating d0{x) to dXui{x) (A parametrizing the associ¬ 
ated zero-energy mode) yields A = f dX = f dx\/2 = 
a/2/. 

Eor the ratio of the determinant r = det' H/| det' in 
Eq. (S3), it is straightforward to generalize the method 
used by McGumber and Halperin [SI] for superconduct¬ 
ing wires owing to the similarity between our eigen¬ 
value problems [Eqs. (S5) and (S7)] and their problems 


[Eqs. (3.12) and (3.13) in the reference]. Applying the 
method of Ref. [SI] with proper boundary conditions [S 8 ] 
yields the exact value of the ratio: r = 4. This value can 
also be numerically verified by evaluating the expression 
with the exact eigenvalues that we are able to obtain with 
the aid of supersymmetry of Eqs. (S5) and (S7) [S9]. 
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